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Abstract 

Effects of quantum statistics are clearly seen in the final states 
of high-energy multiparticle production processes. These effects are 
being widely used to obtain information about the regions where the 
final state hadrons are produced. Here we briefly present and discuss 
the assumptions underlying most of these analyses. 

PACS 25.75. Gz, 13.65,+i 
Bose-Einstein correlations. 



1 Introduction 

Differences between quantum and classical (Boltzmann) statistics show up 
in a variety of systems. Here we will discuss systems, where the energy per 
particle is of the order of 100 MeV or more. In such systems, produced 
in high energy scattering processes, correlations due to quantum statistics 
are clearly seen in the data: Identical bosons seem to attract each other, 
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identical fermions seem to repel each other. The quantitative description of 
these effects is interesting both for their own sake and in order to disentangle 
other more subtle correlations. There is still another motivation for this 
study, however, which is somewhat controversial, but very stimulating. 

Much work has been done on quantum statistics in multiple particle pro- 
duction processes. For reviews see e.g. [T] and jSJ. According to most of this 
work, the study of interparticle correlations due to quantum statistics yields 
valuable information about interaction regions, i.e. about the regions where 
the hadrons are produced. The interaction regions are difficult to study, be- 
cause they are both small and short-lived. Their typical sizes are of the order 
of fermis and also their life-times are of the order of fermis. A fermi in time 
is the time necessary to cross the distance of one fermi at the speed of light 
i.e. about 3 x 10 -24 sek. Very few methods for studying interactions regions 
are available. 

An added attraction is that the predictions from the recognized theories 
for the new data from the heavy ion collider RHIC turned out to be com- 
pletely wrong (references can be traced e.g. from [3]). The name coined 
for this disaster is the RHIC puzzle. Consequently, both the general theo- 
retical framework and the specific phenomenological assumptions have to be 
reanalyzed. In this paper we will present the basic assumptions of the most 
popular models. 

In low energy scattering, i.e. when the center-of-mass kinetic energy of 
the two colliding particles is of the order of 1 GeV, when several particle are 
produced their distribution is roughly spherically symmetric. The first model 
to gain wide popularity, Fermi's model, used a slightly modified microcanon- 
ical distribution. For each multiplicity of particles, each final state with total 
energy and total momentum equal to their initial values was assumed to be 
equally probable. At higher energies spherical symmetry breaks down and 
the model does not work any more. At low energies, however, it seemed good 
until paper jlj got published. 

In this paper the authors studied the opening angles i.e. the angles be- 
tween the momenta for pairs of pions. From Coulomb interactions one would 
qualitatively expect that for pairs of like-sign charged pions the opening an- 
gles would be on the average smaller than for pairs of unlike-sign charged 
pions. Quantitatively, however, one can easily estimate that this effect is 
rather small. The authors expected, nevertheless, a similar effect, because it 
had been predicted that in the tt + 7T~ system there is a resonance, now known 
as the p meson, while for like sign pairs no resonance had been expected. To 
their surprise, the experimental result was just the opposite. The opening 
angles for pairs of like-sign pions tended to be smaller. The result got ex- 
plained in the seminal paper [H] as a result of the Bose-Einstein statistics of 
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the identical pions. For many years the effect had been known as the GGLP 
effect. We will review the GGLP paper in the following section, but now let 
us discuss high energy scattering. 

We will consider central heavy ion (e.g. gold-gold) scattering at center- 
of-mass energies of the order of 100 GeV per pair of colliding nucleons. In 
such collisions, in the center-of-mass system, due to Lorentz contraction, 
both nuclei take the form of thin pancakes. When the two pancakes fly 
through each other there are many nucleon-nucleon interactions, but at this 
energies the directions of flight of the nucleons change little and the pancakes 
survive. When they fly away from each other, many strings are formed and 
stretch connecting the color charges in one nucleon with the opposite color 
charges in the other. These string are a characteristic feature of quantum 
chromodynamics. When two opposite electric charges interact, the well- 
known field of forces extends over all space. In quantum chromodynamics 
the corresponding field is confined to a thin (diameter of the order of one fm) 
tube with the opposite color charges at the two ends. In a high energy central 
heavy ion collision many strings are produced and exist simultaneously. It 
is plausible that they merge and produce a well-defined, roughly cylindrical 
region with the two pancakes at its ends. This region is presumably first 
filled with quarks, antiquarks and gluons. Only later hadrons, mostly pions, 
emerge from it. 

There is a number of questions one would like to ask. What is the size 
and shape of this region? What is its life-time and what is the duration of 
hadronization? Note that the life-time and the duration of hadronization 
are in general different. For instance, the region could live for 8 fm without 
emitting hadrons and then emit all the hadrons within 2 fm. If the content 
of the region can be considered a phase, it would be interesting to know 
its equation of state. This would be valuable information for cosmological 
models of the early stages of the Universe expansion and for models of the 
interiors of neutron stars - perhaps some of them have quark-gluon cores. 
Another set of questions is about the transition of this stuff into hadrons: 
is it a phase transition or a cross-over? If it is a phase transition, is it first 
order or continuous? If it is first order, what is the latent heat? 

2 GGLP or HBT or BE correlations 

Let us start this section with a remark about terminology. Starting from 
the seventies, the original acronym GGLP was being gradually replaced by 
the acronym HBT in honor of Hanbury-Brown and Twiss, who several years 
before the GGLP paper had made a related discussion in astronomy p. 
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They used successfully the Bose-Einstein correlations among photons to mea- 
sure the diameters of distant stars. Now the name GGLP is hardly ever used. 
Kopylov and Podgoretskii explained [Zj how, starting from a more general 
formulation, one can obtain the GGLP case by making a parameter tend to 
zero and the HBT case by making the same parameter tend to infinity. Thus, 
GGLP and HBT are two different limiting cases. In fact, some people replace 
HBT by the more neutral BEC standing for Bose-Einstein correlations. 

In order to facilitate comparison with subsequent developments, we will 
present the GGLP results using a more general notation. GGLP start with 
an input single particle density operator - not to be confused with the actual 
single particle density operator for the system - in the form 

Pi = J d 3 x\x)p(x)( x|. (1) 

This corresponds to particles being produced incoherently from point sources. 
Function p(x) gives the distribution of the sources in space. In this model 
the full information about the size and shape of the interaction region is 
given when function p(x) is known. The corresponding input single particle 
density matrix reads 

P/(p;p') = /p(x)e^ x , (2) 

where q = p — p'. Note that given the density matrix p/(p;p') one can 
unambiguously obtain the distribution of sources p(x) just by inverting the 
Fourier transformation. Unfortunately, this nice feature of the theory will be 
lost, when the theory is made more realistic. 

Interpreting the input density matrix as a density matrix, one would 
obtain the single particle momentum distribution 

<Mp)=Mp;p) = /^p(x) = i. (3) 

This is obviously unrealistic. E.g. very large momenta are forbidden by 
energy conservation. GGLP, however, interpreted Oj(p) as a weight for the 
states allowed by energy and momentum conservation. Thus, the result just 
means that the single particle distributions should be calculated from Fermi's 
model. 

For pairs of identical particles, in general, 

fi 2(Pi, p 2 ) + P/(Pi! Pi)P/(P 2 ; P2)) ( 4 ) 

because the right-hand side does not have the right symmetry with respect 
to the exchange of p : and p 2 . GGLP symmetrized it to get 
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^(pi, p 2 ) = p/(pi; Pi)p/(p 2 ; P2) + p/(Pii p 2 )p/(p 2 ; Pi)- ( 5 ) 



For instance, substituting for the density distribution of sources 



P(x) 



1 




(6) 



one finds 



fi 2(Pi,P 2 ) = 1 + e 



R 2 q 2 



(7) 



Again, as a distribution of momenta this is untenable, but interpreted as a 
weight it nicely reproduces the enhancement at small momentum differences 
seen in the data. 

The GGLP paper has been very influential. It still is the most quoted 
paper in the field. It has, however, its weak points. On the technical side, 
it is easy to calculate the weights, but then the integration over momentum 
space is needed. In the usual case, when neither the nonrelativistic approxi- 
mation, nor the ultrarelativistic approximation (all masses tending to zero) 
is justified, this is a cumbersome task. In practice, as far as I know, no one 
got with this approach beyond three identical particles, while in heavy ion 
collisions, at high energies hundreds of identical particles are being produced. 
Moreover, in order to calculate the momentum integral, one needs the exact 
numbers of all the other particles produced and this is usually not available. 
A reasonable approximate solution of this problem has been found and will 
be presented in the next section. 

Also the physics behind the model is doubtful. Time does not appear 
explicitly. As easily checked, this corresponds to the assumption that all the 
identical particles are produced instantly and simultaneously at some time 
to- This is a most unlikely scenario. Just reflect at the question, in which 
reference frame this assumption should be satisfied. Moreover, classically one 
would like to find the probability distribution for a particle being produced 
at point x with momentum p. In order to get consistency with quantum me- 
chanics one has to compromise. In the GGLP model the particle is produced 
at point x and, in agreement with quantum mechanics, its momentum can 
be anything with equal probability. This is then modified by rejecting the 
states forbidden by energy and/or momentum conservation. It is much more 
likely, however, that the particles are produced as wave packets with some 
finite variances of position and momentum. 
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3 Two particle (reduced) correlation functions 



Consider a two-particle correlation function 

^ 2 (Pi,P 2 ) 



C 2 (Pi,P 2 ) 



^2bckg{Pi, P2) 

Here Q2 is the experimental two-body momentum distribution and &2bckg is 
the corresponding distribution with the Bose-Einstein correlations switched 
off. The latter cannot be obtained from the data without further assump- 
tions, but the experimental groups have various methods for getting reason- 
able approximations to it. It would be more in agreement with the termi- 
nology used in statistics to put into the denominator the product of single 
particle distributions instead, but the correlation function defined here is 
more convenient to study separately the correlations due to quantum statis- 
tics. Assuming that this correlation function can be calculated from the 
GGLP approach without phase space integrations, one finds 

n 2 (Pi, ft) = JWPi. Pa) f 1 + MPi'PJI 2 ) . (9) 

V p(Pi>Pi)p(p 2 >p 2 )y 

From now on we skip the subscript I, though p, strictly speaking, is not 
quite the single particle density matrix of the system. Formula Q, which 
has been proposed by Kopylov [B], is plausible and can be easily used for 
comparison of models with experiment. One measures Q2 and divides it by 
the estimated &2bckg to obtain the experimental C 2 . This is compared with 
the C 2 calculated from the model. The procedure is approximate, but it 
avoids the integration over momentum space. Thus, the information how 
many and what kinds of particles have been produced is not needed. 

For instance, for the Gaussian input density © the calculated correlation 
function is 



c 2 ( Pl ,p 2 ) = i + e - fl2( i 2 . (10) 

Comparing it with the data one obtains the root-mean-square radius of the 
interaction region R. Here by assumption this region is spherically symmet- 
ric. A natural generalization [H], [ID] is to replace the product R 2 q 2 in the 
exponent by q\R\ + q 2 R 2 + q 2 R 2 , where q^, q Q and q s are, respectively, the 
components of q along the beam axis and along some two axes perpendic- 
ular to it. Since the maximum length of the strings, achieved just before 
they break, increases with energy, one would expect for central collisions at 
high energies Rl 3> R ~ R s - Experimentally it is observed that Rl is 
comparable with R Q and R s . This has been explained as follows [TH]. The 
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correlations due to Bose-Einstein statistics are visible only for values of |q| 
of the order of or less. Therefore, they can be used only to measure the 
size of the region where particles with similar momenta are produced. This 
region, named by Sinyukov the homogeneity region, is in general smaller than 
the total interaction region. 

The correction of the physical assumptions is much harder and will be 
discussed in the following sections. Essentially there are two strategies: the 
method of wave packets, with its close relative the method of covariant cur- 
rents, and the methods inspired by the concept of the Wigner function. This 
leaves aside the approaches where a complete model is proposed, which can 
be used to calculate anything. In particular it may be used to calculate the 
correlations among the identical particles at small |q| whether or not they 
have something to do with some interaction region. For an example in this 
category cf. e.g. [TTj . 

4 Wave packets and covariant currents 

Let us replace the GGLP single particle input density operator by 

p = J d A x s J d i p s \ip XsPs )p(x S} p s ){ilj XaPs \ } (11) 
which corresponds to the single particle input density matrix 

p(p;p') = / d 4 x s J d 4 p s ilj XsPs (p)p(x s ,p s )ilj XsP3 (p'). (12) 

Here the time dependence of the density matrix is not explicitly written. 
Function p(x s ,p s ) is the distribution of the space-time four- vectors and the 
energy-momentum four-vectors defining the sources. The distribution of par- 
ticles, however, is consistent with quantum mechanics because of the wave 
functions ip XsPs - There are various ways of using this scenario. 

Kopylov and Podgoretskii, who introduced it [T2*] . assumed that the 
sources differ only by their positions in space-time. Thus 

</Wp) = e^0(p), (13) 

where the fourth component of p, necessary to calculate the product px s , 
is given by p = y/p 2 + m 2 and thus, is not an independent variable. The 
density matrix is 

P(p; PO = 0(P)0V) / d 4 x s J d A p sP (x s ,p s y^ (14) 

and 
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fi 1 (p) = l0(p)l 2 . 



(15) 



This is the first success. The model can reproduce perfectly any single par- 
ticle momentum distribution. On the other hand 

C 2 ( Pl ,p 2 ) = l + |(e^>| 2 , (16) 

where the averaging is over the distribution of sources p(x s ,p s ). Here func- 
tions 0(p) cancel. One way of calculating the distribution p(p s ,x s ) is to 
start with some initial distribution and then to propagate it finding from 
some classical equations, e.g. Newton's or Boltzmann's, the functions x s (t) 
and p s (t). 

Another variant, known as the method of covariant currents is to 

put 

VW S (P) = ^ Xs J (—) , (17) 
\m s J 

where m s is the mass of the source usually put equal to the particle mass. 
In the rest frame of its source each current reduces to the same function 
j(po). In this approach the assumption of Kopylov and Podgoretskii that 
each source yields particles with the same momentum distribution in some 
common frame, e.g. in the center-of-mass frame, is replaced by the more 
plausible assumption that the momentum distribution for particles from any 
source looks the same in the rest frame of this source. 



5 Wigner functions and their generalizations 

In statistical physics there is a well-known method of including simultane- 
ously positions and momenta. One uses the Wigner function related to the 
density matrix in the momentum representation by the formula |14| . [TK] 

^ (X ' K) = / Wf P{K + l q ' K ~ (18) 

where 

K= l -{p + p')- X = ^(x + x'). (19) 

Note that, for further use, K and X are defined as four- vectors, but in 
the definition of the Wigner function only three-vectors and a fixed time 
argument (not written explicitly) appear. The Fourier transformation can be 
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inverted, so that there is a one to one relation between the density matrix and 
the corresponding Wigner function. The introduction of the Wigner function 
solves the problem when the production of all the particles is simultaneous at 
some time t. When the particles are produced in a time interval £2], but so 
that the particles produced at different instants of time do not interfere, one 
can overage over time, with suitable weights, both sides of the equation. The 
time averaged density matrix, which one could hope to measure, is related 
to the time averaged Wigner function, but the important information about 
the time distribution of the production process is lost. When the particles 
produced at different instants of time interfere, even more information is 
lost. It is possible to define objects related by Fourier transformations to the 
various components contributing to the density matrix, but from the point 
of view of interpretation they are very different from Wigner functions |16j . 
Formally, one can write ^1] 

p(p;p') = J d 4 Xe iqX S(X,K), (20) 

where function S(X, K) is known as the emission function. One of the dif- 
ficulties is that the four-dimensional Fourier transformations cannot be in- 
verted. The reason is that for fixed K, from Kq = (p 2 — p' 2 )/2 = one 
finds 

<?o = -j^, (21) 

while in order to invert the Fourier transformation one needs p(p;p'), at 
given K, for all values of q and go- 111 fact, this difficulty is general and 
appears as soon as we introduce time-dependent sources. Thus, there is an 
infinity of emission functions S(X, K) which yield the same density matrix 

p(p;p')- 

For instance one could put 

S(X, K) = 5(X - t)e- iqoXo W(X, K; t). (22) 

This formula is correct in the sense that, as easily seen using the inverse 
of transformation (JTSJl . it yields the correct density matrix. It is, however, 
completely useless, because guessing this Wigner function is just as hard as 
guessing the density matrix in the momentum representation. Actually, this 
formula has a simple physical interpretation. Let us choose a moment of 
time to 5 when all the final hadrons are already present and do not interact 
any more. From this time on the density matrix for each hadron is 

p(p,p';t) = e-^ t -^p(p,p';t ). (23) 



9 



Thus, it is enough to know p(p, p'; to) to predict all that happens later. What 
had happened earlier is irrelevant for this prediction. One can just as well 
assume that all the hadrons got created at t = to- Physically this is a stupid 
assumption, but formally it is good enough to predict the states of the system 
at later times. 

One could ask what is the "reasonable" emission function corresponding 
to the realistic picture of particle emission. A formula for this function has 
been proposed by Shuryak [T2j. Let us consider first the source labelled i. 
If it created particles in the pure state Aj(p), the density matrix for the 
particles produced from this source would be Ai(p)A*(p'). Averaging over 
the parameters of the source, one obtains the density matrix for the particles 
produced by source i without the assumption that the state is pure. Summing 
over the sources one obtains the overall single particle density matrix 

p(p,p') = £(^(p)4V)>- (24) 

i 

Expressing the amplitudes Aj(p) in terms of their Fourier transforms Ji(x) 
one can write for source i 

(A(p)A*(p')> = jd'xj d*Ye^ KY (UX + l -Y)J*(X - l -Y)), (25) 
where Y = x — x'. Summing over i and introducing the notation 

(J(X + \y) J*(X - \Y)) = J2{Ji(X + ~Y)J*(X - \Y)) (26) 

i 

one finds 

p(p,p') = jd'xj d 4 Ye^ KY (J(X + \y)J\X - l -Y)) (27) 

Comparing this formula with the (ambiguous!) definition of the emission 
function (|20[) one finds 

S(X, Y) = J d A Ye lKY (J(X + \Y)J*(X ~\y)). (28) 

The strategy is to guess the emission function, with some free parameters, 
using all the available information about the sources. Then, one calculates 
the density matrix and the correlation function C2 in order to fix the param- 
eters by comparison with experiment. The weak point is that, since there 
is an infinite variety of different emission functions which all give the same 
density matrix, a good fit to the data does not necessarily mean that the 
models being used, and consequently the parameters obtained, make sense. 
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6 Conclusions 



The single particle momentum distribution gives the diagonal element of 
the single particle density matrix. Under some simplifying assumptions the 
correlations due to quantum statistics give information about the out of 
diagonal elements of this matrix. Even if we could find from experiment 
the complete single particle density matrix in the momentum representation, 
which is not quite the case, this is not enough to find unambiguously the 
space-time distribution of the sources. On the other hand, given a good 
model, it is possible to fit the experimental data related to the density matrix 
and thus find some parameters of the model. 

The classical full description: for each particle find its momentum and the 
point in space-time where the particle has been produced, cannot be achieved 
respecting the rules of quantum mechanics. Two popular compromises are: 
use the positions and momenta of the sources instead of the positions and 
momenta of the particles or use X and K instead of x and p. Given the 
general framework one needs some phenomeno logical assumptions to fill it. 
An example, admittedly not a very realistic one, is to assume that all the 
particles are produced simultaneously from a Gaussian distribution of point 
sources. In practice, of course, much more sophisticated phenomenology is 
being used. 

For quantitative work one must include a number of complicating features, 
which have been ignored here for the sake of clarity and of economy of time. 
We give below a partial list. More information can be found in the reviews 

n. 

• Many body symmetrization. 

• Final state interactions. 

• Production of resonances. 

• Momentum-position correlations. 

• Evolution of the interaction region before and during hadronization. 

References 

[1] U.A. Wiedemann and U. Heinz, Phys. Rep. 319(1999)145. 

[2] T. Csorgo, Heavy Ion Physics, 15(200)1. 

[3] K. Zalewski, Nucl. Phys.B (Proc. Suppl.) 133(2004)31. 



11 



[4] G. Goldhaber et al., Phys. Rev. Lett. 3(1959)181. 

[5] G. Goldhaber, S. Goldhaber, W. Lee and A. Pais, Phys. Rev. 
120(1960)300. 

[6] R. Hanbury-Brown and R.Q. Twiss, Nature 178(1956)1046. 

[7] G.I. Kopylov and M.I. Podgoretskii, Zhur. Exp. Theor. Phys. 
69(1975)414. 

[8] G.I Kopylov, Phys. Lett. B50(1974)472. 

[9] P. Grassberg, Nucl. Phys. B 120 (1977) 231. 
[10] G.F. Bertsch, Nucl. Phys. A498( 1989) 173c. 
[11] B. Andersson and W. Hofmann, Phys. Lett. B169(1986)364. 
[12] G.I. Kopylov and M.I. Podgoretskii, Yad. Phys. 19(1974)434. 
[13] M. Gyulassy, K. Kaufmann and L. Wilson, Phys. Rev. C20( 1979)2267. 
[14] S. Pratt, Phys. Rev. Lett. 53(1984)1219. 
[15] A. Bialas and A. Krzywicki, Phys. Lett. B354( 1995) 1345. 
[16] K. Zalewski, Acta Phys. Pol. B34(2003)3379. 
[17] E. Shuryak, Phys. Lett. 44B(1973)387. 
[18] M.G. Bowler, Z. Phys. C29(1985)617. 



12 



